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Abstract 

We demonstrate the relevance of complex Gaussian stochastic processes to the 
stochastic state vector description of non-Markovian open quantum systems. 
These processes express the general Feynman- Vernon path integral propagator 
for open quantum systems as the classical ensemble average over stochastic 
pure state propagators in a natural way. They are the coloured generalization 
of complex Wiener processes in quantum state diffusion stochastic Schrodinger 
equations. 
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I. INTRODUCTION 



The reduced density operator of a quantum subsystem is obtained from the total density 
operator by tracing over the environmental degrees of freedom. Feynman and Vernon |l ] 
derive the propagator J(xj, x'f, t; x , x' Q , 0) of the reduced density matrix p(x, x', t) in terms 
of a double path integral [see also Feynman and Hibbs 0] 

J(x f ,x' f ,t-x ,x' ,0) = r f,t V[x] f Xpt V[x']exp{Us[x]-S[x'])\ T[x,x\ (1) 
1 Jxofl Jx' ,o In J 

where S[x] is the classical action functional of the subsystem alone. The influence functional 
!F[x, x'} combines the effects of the environmental initial state, its Hamiltonian and the 
interaction Hamiltonian between subsystem and environment, on the subsystem. 

It is also shown in [|T],|2| that the general influence functional which is at most quadratic 
in the coordinates must be of the form 

T\x, x'\ = exp j— J dr J da [x T — x' T ][a(r, a)x a — a*(r, c)x^.]| , (2) 

with a positive, Hermitian kernel 

a(r,a) = a*(a,r). (3) 



A. Real kernels 

Feynman and Vernon emphasize that if the kernel is not only Hermitian but real, 
the influence functional can be obtained from a real Gaussian stochastic process F(t) (a 
fluctuating force), with statistical properties 

(F(r)> = (4) 
(F(r)F(a)) = a(r,a) = a*(r,a). 

Here, and throughout the paper, (...) denotes the classical ensemble average over the 
stochastic processes. Using the general formula 
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(exp f dr[f(r)F(r))) = exp ~ f dr f da[f(r)a(r, a)f(a)} (5) 
Jo 2 Jo Jo 

for arbitrary functions /(r), the propagator J for the density operator can be stochastically 
decoupled into stochastic pure-state propagators Gp, 

J(x f ,x' f ,t;x ,x' ,0) = (G F (x f ,t;x ,0)G* F (x' f ,t;x' ,0)), (6) 

with 

Gf(x/, t; xq, 0) = / X>[x] exp j — S[x] + i [ dTX T F{r) \ . (7) 
Jx ,o in Jo J 

The total action functional in the exponent of the path integral propagator (0) now represents 
the additional influence of the stochastic force hF(r) on the subsystem. Thus, real kernels 
are equivalent to ordinary unitary, but stochastic, quantum dynamics. 



B. Complex kernels 

Genuine environments, however, not only induce fluctuations in the subsystem but also 
dissipation. These are represented by complex kernels a(r, a) and can therefore not be 
simulated by a stochastic potential. As an example, Feynman and Vernon derive the 
influence functional (§) analytically for the case of a linear coupling to a heat bath of har- 
monic oscillators with temperature T. Caldeira and Leggett further elaborate this approach 
in [see also Grabert || and Weiss 0], resulting in the complex quantum Brownian motion 
kernel 

a(r, er) = / du uj I coth ( — — J cos (oj(t — a)) — i sin (lu(t — a)) 1 , (8) 



nh Jo [ \2kT 

where £1 is a bath cut-off frequency, m the mass of the particle and 7 the damping rate. 

Remarkably, it has been shown only recently by Diosi || that even in the general case of a 
complex kernel like @, the Feynman- Vernon propagator ([]]) allows a stochastic decoupling 
similar to (^). His result is based on the tricky construction of a real Gaussian process, 
whose correlation function has been given implicitly in terms of a. From the point of view 



of applications, however, the use of these processes appears rather difficult. The deeper 
reason behind this construction comes from relativistic measurement theory 0. 

It is the aim of this paper to present an alternative and much simpler stochastic de- 
coupling of the general Feynman- Vernon influence functional, based on complex Gaussian 
stochastic processes. In their white noise version, these processes have been introduced from 
symmetry considerations in the quantum state diffusion (QSD) stochastic Schrodinger equa- 
tion > describing Markovian open quantum systems. They also appear in measurement 



theories, particularly in cases where the apparatus is represented by a Bosonic reservoir [see 
|13| for more references]. 

Markovian stochastic state vector methods have proven indispensable for many appli- 



cations, particularly in quantum optics []14|] . Our result represents a first step towards an 
applicable non-Markovian stochastic state vector theory, required for instance in solid state 
theory §§. 



In a recent related work Kleinert and Shabanov |TjJ, derive operator quantum Langevin 
equations in the Heisenberg picture corresponding to the propagator (JJ). In this paper, 
however, we stick to path integrals and state vectors in the Schrodinger picture. We review 
basic properties of complex Gaussian processes in Sect. 2. In Sect. 3 we show how they 
enable the stochastic decoupling of the general Feynman- Vernon propagator, resulting in 
linear non-Markovian quantum state diffusion. We conclude with a short discussion and 
further comments in the final Sect. 4. 

II. COMPLEX GAUSSIAN STOCHASTIC PROCESSES 

Here we review complex Gaussian processes Z(r) with stochastic properties 

(Z(r)) = 0, (9) 
(Z(t)Z(o-)} = and 
(Z(r)Z*(a)) = 1 (r,a). 
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Such processes Z[f) can only be constructed if the complex correlation 7 is Hermitian and 
positive, which is automatically fulfilled by the quantum environments we are interested in. 
The processes Z(r) with properties (|9|) are the natural coloured generalization of complex 
Wiener processes £(r) with corresponding Ito increments d£ with properties 

(df) 2 = and I^CI 2 = dt, (10) 

as they arise from symmetry considerations in the quantum state diffusion theory of Marko- 
vian open quantum systems Hq-p] . Complex processes with properties @ are also common 
in quantum measurement theories Jl3| . 



Writing Zir) = X(j) + iY{r) we find that conditions (§) are fulfilled for 

(X(r)X(a)) = (Y(r)Y(a)) = ^Re { 7 (r, a)} and (11) 
(X(r)Y(a)} = -(Y(r)X(a)} = -~Im { 7 (r, a)} . 

We see that the crucial advantage of complex processes is to allow a non- vanishing correlation 
between their real and imaginary part, regarding the real processes X(r) and Y(r) as one 
joint real Gaussian process (X(t),Y(t)). This construction leads to the imaginary part of 
the correlation function 7(1", a) in ([]). 

The relevant formula for complex processes Z{r) that replaces formula (|H) for real pro- 
cesses is 

(exp f dr[f{r)Z{r) +g{r)Z*{r)]) = exp f dr f da [/(t) 7 (t, a)g(a)\ , (12) 
Jo Jo Jo 

valid for arbitrary functions /(r) and g{r). Notice that this implies 

(exp [ t dr[f(r)Z(r)}} = l (13) 

J 

in contrast to equation (|]) for real processes. 

If one wants to generate such processes Z{j) numerically, one can use the construction 

Z(t)= J dr 7 i(*,r)e(r), (14) 



where £(r) is an easily generated white complex process and 

dr 7i(£,t)7i(t, s) = j(t, s). (15) 



III. STOCHASTIC DECOUPLING OF THE FEYNMAN- VERNON INFLUENCE 

FUNCTIONAL 

We show how the processes Z(r) naturally lead to a stochastic decoupling of the 
Feynman- Vernon influence functional. We choose processes with the complex conjugate 
of the Feynman- Vernon kernel from (||D as correlation function, 

7 (r,a) = a*(r,a). (16) 

According to formula flT2|) we find 



(exp / dr [x t Z(t) + x' T Z*{r)}) = exp / dr I da [x T a*(r, a)x' a ] (17) 
Jo Jo Jo 

= exp dr da [x t ol*(t, cr)x' + x' a(r, cr)xJ , 
Jo Jo 

where we used (|3]) to get the second line. This last expression is just the part of the influence 
functional (H) that couples x and x' . Notice, however, that they are decoupled on the left- 
hand side of equation (|TT|). We conclude that we can express the propagator of the density 
matrix in the decoupled form 

J(x f ,x' f ,t;x ,x' ,0) = (G z {xf,t;x ,0)G* z {x' f ,t;x' ,0)}, (18) 

with the stochastic path integral propagator for state vectors, 

Gz(xt, t; x , 0) = / V\x] exp \ — Six] + / drx T Z(r) — dr da [x t oc(t, a)x rT ] \ . (19) 
Jxo,o In Jo Jo Jo J 

Thus, the Feynman- Vernon path integral propagator ([!]) for the density operator is equivalent 
to the ensemble of pure state propagators (|19"D. This is the main result of this paper. The 
non-local action functional in ([T^) reflects the non-Markovian nature of the problem. 
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Result fll"8|) with (|19|) allows to describe the non-Markovian dynamics of the subsystem 
in terms of an ensemble of stochastic state vectors 

\Mt)) = Gz(t ] 0M ). (20) 

If we assume an initial pure state 

Po = P^ , (21) 
where we use the notation = \ifj)(4>\ for pure state projectors, we recover the density 



operator at time t by taking the ensemble average according to ([181) , 



P(t) = (iW- (22) 



In Markovian subsystem dynamics @-ll|, the time evolution of the stochastic state 



vectors is given by a stochastic Schrodinger equation. In fact, in the white noise case 

a(t,s) = nS(t- s), (23) 
the path integral propagator fliTS| ) becomes local in time, 

Gt(x f ,t]x o ,0)= r' t V[x}exp{^S{x] + ^ [' dC(r)x T -^ f' drxA , (24) 

Jxofl I II JO 2 JO ) 

with a delta-correlated normalized complex process f (t) = «-3Z(t). This is the propagator 
of the linear quantum state diffusion stochastic Schrodinger equation 

\di/>) = (-~H - ^x 2 ^j \tf))dt + y/^x\ip)d£ (25) 

with complex Ito increments ([101) . The path integral theory of general linear quantum state 
diffusion equations was developed in |12|], where a rigorous definition of stochastic path 
integrals like (|24] ) is given. 

As the non-Markovian generalization of ([24]), the propagator ( |I~9"D represents linear non- 
Markovian quantum state diffusion. 
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IV. CONCLUSIONS 



We use complex Gaussian stochastic processes to find a non-Markovian quantum state 
diffusion theory which is equivalent to the Feynman- Vernon density matrix formulation. Our 
result offers a helpful tool, since state vectors are simpler than density operators. Such a 
reduction in complexity is essential to tackle many realistic problems and is most significant 
numerically, as is well recognized in the Markovian case, for instance in quantum optics 
|nj. To be truly helpful, one must overcome the difficulties arising from the fact that the 
stochastic state vectors \ipz) of (120) are not normalized. Also an effective algorithm to 



propagate state vectors with the non-local path integral (|T^) has to be established. 

From a quantum foundational point of view the question arises in what sense the stochas- 
tic state vectors \ipz) can be related to the non-Markovian dynamics of individual open 
quantum systems. 

In the well established Markovian case, the use of white complex noise arose from sym- 
metry considerations in the space of the environment operators of the corresponding master 



equation H^TTf. For a single, Hermitian environment operator like x in (|25|), however, this 



symmetry argument does not apply. In this paper we have shown that even in this case, 
complex Gaussian processes appear indispensable as soon as non-Markovian features are 
taken into account. These independent arguments underline their relevance to a stochastic 
description of general open quantum systems. 
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